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ABSTRACT 

We construct updated sol ar models with differ ent sets of solar abundances, including the 
most recent determinations by I Asplund et ah I (|2009l) . The latter work predicts a larger (~ 10%) 



solar metallicity compared to previous measurements by the sa me authors but significantly lower 
(~ 25%) than the recommended value from a decade ago by iGrevesse fc Sauval ( 1998f ). We 



compare the results of our models with determinations of the solar structure inferred through 
helioseismology measurements. The model that uses the most recent solar abundance determi- 
nations predicts the base of the solar convective envelope to be located at _Rcz = 0.724Ro and a 
surface helium mass fraction of ysurf = 0.231. These results are in conflict with helioseismology 
data (i?cz = 0.713 ± O.OOIR0 and Fsurf = 0.2485 ± 0.0035) at 5-cr and ll-a levels respectively. 
Using the new solar abundances, we calculate the magnitude by which radiative opacities should 
be modified in order to restore agreement with helioseismology. We find that a maximum change 
of ^ 15% at the base of the convective zone is required with a smooth decrease towards the core, 
where the change needed is ~ 5%. The required change at the base of the convective envelope is 
about half the value estimated previously. We also present the solar neutrino fluxes predicted by 
the new models. The most important changes brought about by the new solar abundances are 
the increase by ^ 10% in the predicted ^^N and ^^O fluxes that arise mostly due to the increase 
in the C and N abundances in the newly determined solar composition. 

Subject headings: Sun: helioseismology - Sun: interior - Sun: abundances - neutrinos 



Introduction 



agreement between standard solar model (SSM) 
predictions and helio seismology determinations of 
the solar structure (IChristensen-Dalsgaard et al.l 



Since the solar surface heavy-e lement content 
has b een revised downwards by lAsplund et al. 

AGS05), from (Z/X)© = 0.0229 (jGrevesse fc Sauval ^.Qj^^ g3(j ( Basu fc Antia 



(1200,'i 

Il998l GS98) to (Z/X)© = 0.0165, the excellent 



9961 : iBahcall et al.ll200lh has b e en seriously com 



2004 iBahcaU et al.l 120051) . This discrepancy be- 



2004 iMontalban et al 
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tween models and helioseismic inferences has 
prompted a number of authors to revise the 
physical inputs of S S Ms Qv lontalban ct al. 200; 



Basu fc Antial |2004 iBahcall et al 



2006 



Christensen-Dalsgaard et al 



I2OO5I : iGuzf 
20091 ). and to 



question the revision of the solar abundances, par- 
ticularly those of C, N, O, Ne and Ar whose frac- 
tional abundances cannot be determined from 



Bahcall et al. 


2005; 


Delahave & Pinsonneaultl2006t 


Basu & Antia 


2008 


). 



Very recently, lAsplund et al. ( 20091 hereafter 



AGSS09) have done a complete revision of the 
solar photospheric abundances for nearly all ele- 
ments, including a new 3D hydrodynamical solar 
atmosphere model with improved radiative trans- 
fer and opacities (^Trampedac h et al lliooi). The 
predictions from this 3D model have been shown 
to agree remarkably well with various observa- 
tional constraints, including the atmospheric ther- 
mal structure as judged from continuum center- 
to-limb variat i on and detailed line profile shapes 
( Pereira et al. 2009al |bh. The newly determined 
solar abundances lead to {Z/X)q = 0.0178, higher 
than the AGS05 value but still well below older 
determinations, e.g. GS98. 

In this Letter we present a series of new SSM 
calculations using solar compositions from GS98, 
AGS05, and the newly determined solar abun- 
dances by AGSS09. All the models incorporate 
new refinements in the input physics, so that 
models presented here with the older compositions 
(GS98 and AGS05) represent updat ed versions of 



previ ous solar model calculations (jBahcall et al 
120051 ). For each model, we compare our results 
with helioseismological determinations of solar 
properties and also give the predicted solar neu- 
trino fluxes. Additionally, and motivated by the 
new AGSS09 composition, we determine the fac- 
tor by which radiative opacities in the solar inte- 
rior should be increased to solve the solar abun- 
dance problem following the scheme presented by 
Christensen-Dalsgaard et al. I llooi). 



2. Calculations and Results 

Solar models in this work have been com- 
puted with a modifie d version of GARSTEC 
(| Weiss fc Schlatt]||2008l) that uses the nuclear en- 



ergy generation routine exportenergy.^. Element 
diffusio n in the sola r inter ior is included accord- 
ing to Thoul et"al] ( 1994h . Radiative opacities 
are from the Opacity Project, co mplemented at 



low te mperatures with those from [Ferguson et al 
(|2005l ). Specific sets of opacities have been com- 



puted for each of the solar compositions used in 
this paper (see below). With respe c t to p revi- 



ous works, e.g. iBahcall et all (|2005L I2OO6I ). the 



changes in the input physics are: a revised version 
of the OPAL equation of stateH (EOS) that cor- 
rects errors in the 2 002 O PAL EOS tables from 
Rogers fc Navfonov ( 2002f ) (our previous choice), 



and updated values of tw o important nuclear as - 
trophysica l factors, S^zi (ICostantini et al.l l2008f ) 



and Si 14 ( Marta et al. 20081 ). the latest determi- 



nations by the LUNA experiment. 

We have computed solar models using three dif- 
ferent basic solar abundances. Two models em- 
ploy previous solar abundance compilations (GS98 
and AGS05) and show small differences with re- 
spect to models with the same abundance s pre- 
sented elsewhere, e.g. IBahcall et al.l (|2005t ). The 
changes originate from the use of the updated EOS 
and cross-sections for nuclear reactions mentioned 
above. A third SSM has been computed adopting 
the new solar composition determined by AGSS09. 
The most important results in this work are re- 
lated to this model. The choice of the abundance 
scale (meteoritic or photospheric) deserves a short 
discussion. While AGSS09 find the average differ- 
ence between photospheric and meteoritic abun- 
dances to be 0.00 ± 0.04 dex, a few elements rel- 
evant to detailed solar modeling show compara- 
ble or slightly larger deviations. This is the case 
for Mg, Ca, and Fe for which differences between 
the two scales are 0.07, 0.05, and 0.05 dex respec- 
tively, photospheric values being larger. Given the 
historical robustness and higher accuracy of me- 
teoritic determinations of abundance ratios, and 
the present excellent overall agreement with pho- 
tospheric abundances, we adopt for AGSS09 me- 



^ Publicly available at http:/ /www. sns.ias.edu/~jnb' 
2http://adg.llnl.gov/Researcli/OPAL/EOS.2005/. OPAL 
EOS us es relative m etal abundances for C, N, O, and Ne 
close to lGrevesse &: N oels (1993) and abundances of heav- 
ier elements are added to Ne. Differences in relative abun- 
dances of metals with more recent solar abundance com- 
pilations have negligible influence in the global properties 
of sol ar models, provided the correct overall met allicity is 
used llGong et alj|200ll : iRoeers fc Navfonovll2002t) . 
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teoritic abundances for refractory elements as the 
standard choice for our solar models. We also in- 
vestigate, however, the use of adopting the photo- 
spheric values throughout. This is also consistent 
with the adoption of meteorit ic scales in previous 



work s on solar models (e.g. Bahcall et al. 20051 



20061 ) . The abundances for the different solar com- 



positions used in this paper are given in Table [l] 
in particular for those elements entering the cal- 
culation of radiative opacities. 

Models have been evolved from the pre-main se- 
quence to the present-day solar system age, tq = 
4.57 Gyid. The main characteristics of the mod- 
els are listed in Table [2] The second and third 
columns give the present-day heavy elements to 
hydrogen mass ratio and the surface metallicity. 
From the fourth to the seventh columns we present 
quantities directly related to hclioseismology: sur- 
face helium mass fraction ygurf , depth of the con- 
vective zone i?cz /R-o , and the average rms differ- 
ence between model and solar sound speed and 
density profiles, {5c/ c) and (Sp/p) respectively. 
Columns eighth and ninth give the central helium 
mass fraction Yc and metallicity Zc at tq . Finally, 
the last three columns give the initial composition 
of the models and the mixing length parameter. 

Results for the GS98 and AGS 05 models are 
very s imilar to those presented in iBahcall et al. 
(|2005l ): the improved EOS leads to changes 



in the sound speed and density profiles about 
one order of magnitude smaller than the dif- 
ferences between solar models and helioseismic 
inferences, while changes in nuclear cross sec- 
tions only affect predictions for neutrino fluxes 
that are discussed in detail below. Helioseis- 
mically derived values for Rcz/K(7) and Fsurf 
and are 0.713 ±0.001 (iBasu fc Antial Il997l) and 



0.2485 ± 0.0035 (|Basu fc Antiall2004l ) respectivelv. 
The GS98 model predicts Rcz/Rq with the right 
value compared to hclioseismology, while the 
AGS05 model shows a 15— cr discrepancy. For 
Ysuri the situation is analogous: GS98 value dif- 
fers from the hclioseismology determination by 
1.8-a while for AGS05 the discrepancy is 5.5-a. 

For the SSM that adopts the newly deter- 
mined solar composition AGSS09, i?cz/R-o and 



l^urf show some improvements with respect to the 
AGS05 model, but still far away from the hclio- 
seismology values by ll-cr and 5-a respectively. 
The slight increase in oxygen abundance (0.03 dex) 
and the larger change in neon (0.09 dex) contribute 
to enhance the opacity below the convective zone 
(CZ), deepening its inner boundary and decreas- 
ing the mismatch with the solar sound speed at 
the same time. This is illustrated in the top panel 
of Figure [T] where the relative difference in sound 
speed is shown for the models considered in this 
work. The peak in the profile of the sound speed 
difference, right below the CZ, is ^ 1% for the 
AGSS09 composition, an improvement with re- 
spect to AGS05 but still significantly higher than 
that for the GS98 model. The overall agreement 
with the solar sound speed, as derived by inver- 
sions and indicated by (Sc/c), is a factor of four 
times worse for the AGSS09 model than for GS98. 
For the latter model, (Sc/c) is almost unchanged 
compared to results from previous works. How- 
ever, a detailed comparison of the sound speed 
pr ofiles shown in Figu re [T] with those presented 



•^Detailed structure of solar models at tq 
presented in this work can b e found in 
[http: / /www.mpa-garching.mpg.de/~aldos| 



in iBahcall et al.l (|2005l ) unveils some differences, 
more evident below R ~ 0.6 R©. These changes 
are not related to the improved physics adopted 
in the new models, but rather the result of us- 
ing better data for low-degre e (£ < 3) modes t hat 
penetrate the solar core (see iBasu et al.l l2009l for 
details). In the same figure, results for density 
inversions are shown in the bottom panel. Again, 
the AGSS09 composition leads to an improvement 
in the agreement with hclioseismology compared 
to the AGS05 model, but still far from the results 
obtained for the GS98 composition. 

As already mentioned, meteoritic and photo- 
spheric abundances in AGSS09 agree with each 
other very well, but a few elements show differ- 
ences that could have potential impact on details 
in the solar structure. To quantify this assertion, 
we have computed an additional SSM using only 
the photospheric abundances given in AGSS09, for 
which {Z/X)q — 0.0181. The main characteristics 
of this model, identified as AGSS09ph, are given 
in the last entry of Table [21 Compared to the 
model with the meteoritic abundances, AGSS09ph 
performs somewhat better in terms of helioseis- 
mological quantities as inferred from the results 
summarized in the table, with discrepancies with 
the measured depth of the CZ and surface helium 
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abundance of the order of 9-a and A-a respectively. 
The sound speed and density profiles are shown as 
dotted lines in Figure [T] The changes with respect 
to our standard AGSS09 (meteoritic scale) model 
changes are mostly due to the larger Mg and Fe 
photospheric abundances (0.07 and 0.05 dex, re- 
spectively) that enhance the opacity in the radia- 
tive interior; the fractional increase in opacity is 
larger than the fractional increase in the overall 
metallicity (note the largest improvement in the 
sound speed, for instance, occurs at i? w 0.5 R©, 
the region where the contribution of Mg to the 
opacity is largest). Our results show that adop- 
tion of the photospheric scale gives slight improve- 
ments in the solar model predictions. However, 
since the uncertainties in the meteoritic abun- 
dances typically are smaller than the correspond- 
ing ones for the photospheric determinations the 
meteoritic scale is our preferred choice for solar 
abundances (with the exception of the volatile el- 
ements that are depleted in meteorites) . This is re- 
inforced by the historical robustness of meteoritic 
abundance determinations. 

Low-degree helioscismology {£ < 3) can be 
used to derive seismic information on the solar 
core, particularly by usi ng the so-called separa- 



tion ra ti os as described bvlRoxburgh fc Vorontsov 



(|2003l ). IChaplin et all (l2007|) have used the sep- 



aration ratios constructed with very precise fre- 
quencies of low-^ modes measured by the Birming- 
ham Solar-Oscillations Network (BiSON) to con- 
strain properties of the solar core. They showed 
the discrepancy in solar models constructed with 
AGS05 composition extends all the way to the so- 
lar core and it is not just related to deficiencies in 
the modeling of the solar outer layers, in the con- 
vective envelope. Here we compare the observed 
separation ratios to those computed for our solar 
models; results are displayed in Figure [21 As with 
other helioscismology indicators, the GS98 model 
performs much better than the AGS05 model. The 
adoption of the AGSS09 composition in the SSM 
has very small influence in the core structure of the 
model, as it practically overlaps the AGS05 model. 
Results for the AGSS09ph model closely resemble 
those from AGSS09 and, for cl arity, are no t show n 
in Figure [H As discussed in Basu et al] ( 2007 ). 
values of the separation ratios are closely related 
to the quantity 1/r [dc/dr) integrated over the so- 
lar structure. Differences in this quantity between 



models with AGS05 and AGSS09 (both meteoritic 
and photospheric) compositions are very small and 
only present very close to the center [r < 0.05 i?©) 
and thus have a very small impact on the separa- 
tion ratios. 

Solar neutrino fluxes for the models are listed 
i n Table [3l The new astrophysical factor S34 



( Costantini et al. 20081 ) is 7% larger than previous 
determination and it is responsible for the larger 
^Be and ®B fluxes of models GS98 and AGS05 with 
respect to those published in Bahcall et al. ( 20061 ) 
with the same solar compositions. Similarly, the 
somewhat smaller (7.6%) 5*1 _ 14 value recently pub- 



lished by the LUNA collaboration (jMarta et al 



20081 ) leads to a proportional reduction in the N 



and "'^^O fluxes. The increase in metallicity in the 
AGSS09 solar composition, compared to AGS05, 
is only partially reflected in the changes of the 
neutrino fluxes. This is so because elements that 
more strongly affect the solar core temperature 
(Si, S, and Fe) have the same abundance in AGS05 
and the meteoritic AGSS09 scale. The increase 
by 0.09 dex in Ne and by 0.22 dex in Ar are the 
most important changes in abundances influenc- 
ing the core temperature. This is reflected, for 
instance, in the ~ 4% increase in the *B flux from 
AGS05 to AGSS09. and fluxes, that de- 
pend almost linearly in the added abundance of 
carbon and nitrogen show larger changes (of the 
order of - 12 - 14%) due to the 0.03 and 0.05 dex 
larger abundances of these elements in the new 
AGSS09 abundances. Differences in the neutrino 
fluxes between the AGSS09 and the GS98 mod- 
els are of the order of 10% for '^Be, 20% for **B 
and 38% for the added and ^^O fluxes. For 
the sake of completeness we present the neutrino 
fluxes for the AGSS09ph model. The effect of the 
increased iron abundance is readily noticeable par- 
ticularly in the larger ^B flux compared to AGSS09 
model and, to a lesser extent in the CNO and ^Be 
fluxes. The possibilities that current and future 
neutrino experiments offer to constrain the solar 
core metallicity are beyond th e scope of this paper 
and arc discussed elsewhere ("Hax ton fc Serenelli 
2008; Pefia-Garay & ScrcneUi 200^^ 

Qualitatively, the AGSS09 abundances do not 
change the picture that emerged with the previ- 
ous set of solar abundances, AGS05: SSMs con- 
structed with abundances derived from the most 
sophisticated solar atmosphere models and up-to- 
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date atomic data conflict with all helioseismology 
inferences of solar structure. Potential solutions to 
the solar abundance problem analyzed by differ- 
ent authors, none of them successful in restoring 
the agreement between SSM and helioseismology 
measurements, still face the same problems if the 
AGSS09 composition is used instead of AGS05. 
In this r egard, it is particularly informa tive the 
work by IPelahaye k, Pinsonneault ( 20061 ) , which 



show in the i?cz — ^surf plane the direction in 
which model predictions change when modifica- 
tions in their input physics are applied. It is 
not a simple task to find deficits in the mod- 
els that produce simultaneously changes in both 
RcT, and Fsurf in the right direction, except for 
restoring the metallicity to a larger value, com- 
parable to that of GS98. Quantitatively, how- 
ever, the disagreement is less severe with AGSS09 
composition and, motivated by this fact, it is 
worth taking another look at the opacity deficit 
induced by the AGSS09 abundances. We have 
done a similar analysis to t h e one presented by 
Christensen-Dalsgaard et al. (|2009l ) to check by 
how much opacities in the AGSS09 model have 
to be increased to recover the level of agree- 
ment with helioseismology that the GS98 model 
gives. We find that in the central regions the re- 
quired change is ~ 5% (2% with AGSS09ph) very 
close to what IChristensen-Dalsgaard et al.l ( 2009[ ) 
found. The magnitude of the change increases 
smoothly outwards and reaches ^ 15% (12% with 
AGSS09ph) at the bottom of the GZ. This re- 
quirement is smalle r by almost a factor of two 
than that found by iGhristensen-Dalsgaard et aT 



( 2009( ). very likely because they used the S-model 
(which used the higher {Z/X)q = 0.0245 from 
Grevesse fc Noelslll993l) as their reference model 
and one with the AGS05 solar compositi on, while 
we h ave used GS98 and AGSS09 (see ISereneUi 
20101 for a more thorough discussion) . We note the 
required changes are much larger than differences 
found between OP and OPAL opacities in the ra- 
diative solar interio r. It would be interesting, 
as po inted out by IGhristensen-Dalsgaard et al 



(120091 ). ;o find other observable implications that 
changes of 12 — 15% in the radiative opacity at 
temperatures of a couple to a few millon degrees 
would have, independently of those from solar 
models. 



3. Summary 

We have computed new standard solar mod- 
els that incorporate the most up-to-date input 
physics, including the updated OPAL equation 
of state and the most recent determinations 
of the astrophysical factors for the important 
3He('*He,7)^Be and i'*N(p, 7)1^0 reactions. We 
have used three different sets of solar abundances: 
GS98, AGS05 and the newly determined AGSS09. 
We have found that the updated physical inputs 
have very little effect on the properties of solar 
models, both in terms of solar structure and in 
the neutrino fluxes. The most important results 
in this work are those from the model that adopts 
the new set of solar abundances, AGSS09. The 
new abundances are determined by using an im- 
proved solar atmosphere model and atomic data, 
and a more careful selection of spectroscopic lines. 
The resulting solar metallicity is somewhat larger 
than that from AGS05. This is reflected in the he- 
lioseismological properties of the solar model. For 
the AGSS09 model, the sound speed and density 
profiles, the predicted surface helium abundance 
and the depth of the convective zone are still in 
conflict with helioseismology data, although the 
disagreement is less severe than for the AGS05 
model. Still, results are far from the excellent 
match found with the GS98 composition. Finally, 
we have found that if radiative opacities were to be 
modified to restore the agreement between solar 
models (with AGSS09 composition) and helioseis- 
mology, the required changes are « 12 — 15% right 
below the convective zone with a smooth decrease 
towards the central regions, where changes should 
be 2 - 5%. 

AMS thanks Maria Bergemann for insightful 
discussions on the physics of non-LTE effects in 
the formation of spectral lines. SB is partially 
supported by NSF grant ATM-0348837. We thank 
the anonymous referee for useful suggestions and 
comments. 
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Table 1 

Adopted solar chemical compositions. 



loge 



Element 


GS98 


AGS05'^ 


AGSS09'^ 


AGSS09ph'' 


C 


8.52 


8.39 


8.43 


8.43 


N 


7.92 


7.78 


7.83 


7.83 





8.83 


8.66 


8.69 


8.69 


Ne 


8.08 


7.84 


7.93 


7.93 


Na 


6.32 


6.27 


6.27 


6.24 


Mg 


7.58 


7.53 


7.53 


7.60 


Al 


6.49 


6.43 


6.43 


6.45 


Si 


7.56 


7.51 


7.51 


7.51 


S 


7.20 


7.16 


7.15 


7.12 


Ar 


6.40 


6.18 


6.40 


6.40 


Ca 


6.35 


6.29 


6.29 


6.34 


Cr 


5.69 


5.63 


5.64 


5.64 


Mn 


5.53 


5.47 


5.48 


5.43 


Fe 


7.50 


7.45 


7.45 


7.50 


Ni 


6.25 


6.19 


6.20 


6.22 



Note. — Abudances given as logCj = log Ni/NH + 12. 
" The adopted abundances are the recommended solar pho- 
tosphcric abundances for the volatile elements (C, N, O, Ne 
and Ar) and the CI chondritic meteoritic values for the re- 
maining elements. ^ The adopted abundances are the rec- 
ommended solar photospheric abundances throughout. 
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Table 2 

Main characteristics of solar models. 



Model 


{Z/X)smf 




^urf 


Rcz/Rq 


{6c/c) 






Zc 


^ni 


Zini 


aMLT 


GS98 


0.0229 


0.0170 


0.2423 


0.713 


0.0010 


0.011 


0.6330 


0.0201 


0.2721 


0.0187 


2.15 


AGS05 


0.0165 


0.0126 


0.2292 


0.728 


0.0049 


0.048 


0.6195 


0.0149 


0.2593 


0.0139 


2.10 


AGSS09 


0.0178 


0.0134 


0.2314 


0.724 


0.0038 


0.040 


0.6220 


0.0160 


0.2617 


0.0149 


2.09 


AGSS09ph 


0.0181 


0.0136 


0.2349 


0.722 


0.0031 


0.033 


0.6263 


0.0161 


0.2653 


0.0151 


2.12 
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Table 3 
Predicted neutrino fluxes. 



Model pp pep hep '^Be ^B ^^N ^""O ^'^F 



GS98 5.97 1.41 7.91 5.08 5.88 2.82 2.09 5.65 

AGS05 6.04 1.44 8.24 4.54 4.66 1.85 1.29 3.14 

AGSS09 6.03 1.44 8.18 4.64 4.85 2.07 1.47 3.48 

AGSS09ph 6.01 1.43 8.10 4.79 5.22 2.15 1.55 3.70 



Note. — Neutrino fluxes are given in units of 10^"(pp), lO^(^Be), 
108(pep, ^^N, 1^0), 106(«B, i^F) and 103(hep) cm-^ g-i. Direct 
measurement by the Borexino experiment of the '^Be flux gives 
5.18 lb 0.51 X 10^ cm~^ s~^ after 192 days of data (jArpesella et al. 
l2008f) . For *B, and until the SNO collaboration presents a combined 
analysis of the three phases of the experiment, a simple weighe d av- 
erage of the three phases of the SNO experiment (Ahmad et al..i2002l : 
lAharmim et~ani2005l [2008h gives 5.18 ±0.29 x 10^ cm'^ s'^ for this 
flux. 
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Fig. 1. — Relative sound speed Scjc and 
density 5p/p differences, in the sense (Sun - 
Model)/Model, between solar models and helio- 
seismological results. Details on the inversion pro- 
cedure and data used, as well as the reference 



sound speeds and densities are given in lBasu et al 
()2009l) . 



11 




ISOO 2000 2500 3000 HaOO 

Frequency (^Hz) 




lO F. ..... I .... I ......... I . 

1500 2000 2500 3000 3500 

Frequency (//Hz) 



Fig. 2. — Separation ratios. Comparison between 
values determined from BiSON data and the solar 
models presented in this work. 
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